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p^ I Geometric phases are an interesting resource for quantum computation, also in view of 

their robustness against decoherence effects. We study here the effects of the environment 

'^ . on a class of one-qubit holonomic gates that have been recently shown to be characterized 

OO ' by "optimal" working times. We numerically analyze the behavior of these optimal points 

^^ ' and focus on their robustness against noise. 

o' 
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I. INTRODUCTION 
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Quantum algorithms based on geometric phases [l|, |2] are attracting increasing interest in 
^ . quantum computation BQ. The .elated quantum ,a.e. „p.ese„.i„. the unitary fansform. 

C^ ' tions on a register of qubits do not have a dynamical origin: the Hamiltonian depends on time 

through a set of control parameters that change by following suitable closed loops in the asso- 
ciated parameter space; in the adiabatic limit the dynamical contribution to the evolution can 
be factorized and the features of the quantum gate depend only on the topological structure 
of the manifold. 

Geometric quantum computation has been investigated using both Abelian [5| and non- 



Abelian Idl holonomies. There have 



ion traps 7], Josephson junctions 



m 



3een several proposals for their implementation using 



and semiconductors 10]. Since all physical devices 
interact with their environment, one must carefully analyze the onset of decoherence iljj and 
its detrimental effects against realistic implementations of quantum gates and algorithms. In 
particular, the effects of noise for non- Abelian holonomies in open quantum systems have been 



recently investigated in a series of articles 



jyuyu- inU 



we studied a class of 



one-qubit gates implemented on a four-level ("tripod") system J7|, focusing on non-adiabatic 
effects and bringing to light the presence of fidelity revivals, namely an infinite number of 
(optimal) times at which the fidelity reaches unity. 

In this article we shall investigate the behavior of the fidelity at the first of these optimal 
working points and study its robustness against noise effects. The deviations from the ideal 
(noiseless) case will be numerically analyzed as a function of the strength of the noise (the 
coupling of the system with its environment) and a heuristic definition on robustness will be 
introduced. 

This paper is organized as follows. We review the concept of holonomy in Section ^ and 
briefly introduce the specific tripod system [3] in Section IIIIl where we focus on the role of 
non-adiabatic effects. In Section IIVI we outline the main features of the master equation for 
time dependent Hamiltonians: this is numerically solved in Section in order to analyze the 
behavior of the optimal working points in the presence of noise. In Section IVll we and discuss 
the robustness of our gates. 

II. ABELIAN AND NON ABELIAN HOLONOMIES 

We consider a system governed by a non degenerate Hamiltonian that depends on time 
through a set of parameters, adiabatically covering a closed loop in the parameter space. 
Under these conditions, the final state exhibits, in addition to the dynamical phase, also a 
geometric phase, whose structure depends only on the topological properties of the parameter 
manifold [ly]. If the Hamiltonian has some degeneracies, a loop in the parameter space involves 



more complex geometric transformations J20|. We suppose that the family of Hamiltonians 
H{x{t)) {x'^{t) being a set of parameters) is iso-degenerate, i.e. that the dimensions of its 
eigenspaces do not depend on the parameters and the eigenprojections Pmix{t)) (m denoting 
the eigenvalue) have a smooth dependence on t (at least twice continuously differentiable) . 
In particular, this entails the absence of level crossing between different eigenspaces. H{t) 
can be decomposed by using its instantaneous eigenprojections H(t) = J2m^m{t)Pmit), with 
Pm{t) = Z^fc l'^fc(i))('7ifc(i)| and k the degeneracy index. We define the operator i?, that 
transports every eigenprojection from to to i, and its hermitian generator D{t,to), 

R{t,to)Pm{to) = Pm{t)R{t,to), D{t,to) = -iR{t,to)^ -^^R{t,to). (2.1) 

In the adiabatic limit the evolution of the state remains confined in the degenerate eigenspaces 
and the evolution operator U becomes block-diagonal. In the case of cyclic evolution {Pm{t) = 



Pm{to)) 

U{t,to) ~ ^P„.(to)e"*^*o^'-W'^^C/-p„(to), C/,^ = Pexp|-£^-(x)| , (2.2) 

and the geometric evolution is given by a path ordered integral (P in the above formula) of 
the adiabatic connection A^{x) = '^nA'^dx'^, with 

A^{x{t)) = P„^{x{to))RHx{t),x{to))-^R{x{t),x{to))Pm{x{to)). (2.3) 

If the eigenvalues Cm are time-independent and the connection D is piecewise constant (i.e. 
D{t,to) = D{to,tQ) Vs G [i,to]) the evolution operator reduces to the useful expression [ij] 

U{t, to) = e*(*-*o)^(*0'*o) e-iii-io){H(to)+D(to,to))_ (-2^4) 

We will study a large class of gates where the above hypoteses are satisfied and one can exactly 
evaluate the time evolution, including all non-adiabatic effects, by making use of (|2.4|) . 

III. FREE IDEAL EVOLUTION 

We focus on the "tripod" system introduced in l7"| for holonomic quantum computation: see 
Fig. ^a), where three degenerate levels are connected with a fourth one by Rabi oscillations. 
The adiabatic evo 
implementations 

also non-adiabatic effects. At time t = the logical states and 1 are encoded respectively 
in the quantum states |0) and |1), while \a) is an ancilla state used as "buffer" during the 
evolution. The Hamiltonian of the system reads i:f(t) = |e)(r2o(t)(0|-|-r2i(t)(l|-|-r2a(t)(a|)+H.c., 
where il,j{t) represent the time dependent Rabi frequencies of the transitions. The loop in the 
parameter space is obtained by varying 0,j{t) (j = 0,1, a). In our calculations we consider 



ution of this system was analyzed in several articles for different experimental 
3) M, UM- Let us first review the ideal noiseless case, taking into account 



nj{t) £ R,\/i. The eigenvalues of the system are {0, ±y^0o(t)2 + Qiit)"^ + l^a(i)^ = ±^}, 
where is 2-fold degenerate, corresponding to a 2-dimensional (computational) eigenspace, 
and 0, is kept constant. Therefore, the parameter space is the 2-sphere of radius 0,, {fij G 
R| ^ • il? = ri^}, shown in Fig. ^b). Introducing the parametrization 

ill = rj siwd cos ip, il.Q = i} siwd sin 93, il^ = il cost?, (3-1) 

the eigenstates take the form 

\DQ{t)) = cos(^ |0) — sinc^ |1), 

\Di{t)) = COST? sin(^|0) + cost? cos V3|l) — sint?|a), (3.2) 

\D±{t)) = (±|e)+sini?sin(/?|0)+sini?cos99|l)+cos??|a))/(\/2). 
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FIG. 1: (a): Scheme of a tripod system. and 1 are the computational levels, while a is an ancilla 
state used for the intermediate steps of the transformation. The three degenerate levels are connected 
with an upper level e by time dependent Rabi frequencies i^j{t). The yellow wavy line represents the 
noise, that induces additional transitions only between and e. (b): Path in parameter space for the 
realization of a NOT gate. The solid angle spanned during the evolution is 7r/2. 



The computational space Cs = Span{\DQ{t)) , \Di{t))} belongs to the degenerate eigenvalue 0, 
while \D±{t)) are the bright eigenstates belonging to ±0,. One easily shows that for a closed 
loop on the 2-sphere in the computational space the holonomy (|2.2j) reads C/ad = exp(io"yu;), 
where (Ty = — i(|Z)o(*o))(-Di(to)| — |^i(^o))(^o(^o)|) and w is the solid angle enclosed by the loop. 
In particular, if w = tt/2, we have Uj^/2 = exp(20"y it/2) = iay (in the basis {|L'o(io))i l-^i(*o))})) 
that represents a NOT transformation (up to a phase for the state \Dq)). 

Following the discussion of the previous section we discuss the non-adiabatic corrections to 
this system. In order to use Eq. H2.4() . we will consider the loop shown in Fig. E^b), enclosing 
the solid angle it/2; in the adiabatic limit (when J7t ^ oo, r being the total time of the cyclic 
evolution and fi the energy of the bright states) this path yields a NOT gate. The first step 
consists in constructing the operator D from Eq. 1)3. 2() and the definition ()2.1() . One can see 
|14| that, as far as the rate of change of the polar angles is constant in each section of the path, 
D{t, to) is piecewise constant and we can use Eq. 1)2. 4() to evaluate the evolution operator along 
the path shown in Fig. ^^b). 

An interesting feature of the evolution is that it is factorized in three terms. In the adiabatic 
limit it simplifies to 

/ 1 n n \ 

-1 



u^/2inT) = Usinrs) U2{^T2) Uiinn) 



rfi^QO jjod 



u^j.inr) 



e- 
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3+i-rQ 



(3.3) 



\ 
r being the total evolution time needed for covering the loop in the parameter space and 
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FIG. 2: Mean fidelity versus the cyclic time fir (noiseless case). H. is the energy gap between the bright 
and dark states, r is the time needed to cover the loop shown in Fig.^b). The average is performed 
over a set of initial states uniformly distributed on the Bloch sphere. The dot on the first significant 
peak indicates the optimal working point, Fq = 1 and firj* — 18.25. 



Tj = OiT, with ^j Oj = 1. This represents a NOT gate for the degenerate subspace and yields 
(fast oscillating) dynamical phases for the bright states. 

In order to understand how close the evolution operator is to the ideal one, we use the mean 
fidelity 



{F)(nT) = — / d cos M(pF(nT,'d,ip) 



A-K 



dcos i^dipTilasxi^Qr, 'd, {p)a{i^T, -;?, cp)}, 



where 



aad(f^T,^,V9) = U^%inT)a{^,^)U^%inT), 



(3.4) 

(3.5) 
(3.6) 



a{i}, if) = I"!?, </?)('i9, ^\ being the initial state (assumed to be pure). In practice, in our analysis, 
F will always be averaged over a finite set of input states uniformly distributed on the Bloch 
sphere. The mean fidelity is plotted in Fig. EJ as a function of the adiabaticity parameter VLt. 
Clearly, (F) asymptotically approaches unity (with some oscillations), as expected (adiabatic 
limit). Notice that the fidelity is exactly one for some finite values of time, t = t^, that are 
independent of the initial state. In this case the NOT transformation is perfect, even though 
one is far from the adiabatic regime. 



It is possible to show that, when the three arcs in the loop in Fig. ^b) are covered in equal 
times, one obtains 



r* = —Vm^^, kGN*. (3.7) 

The first fidelity revival occurs for k = 1 

Oti* = ^\/l5 = 18.25 (3.8) 

and is indicated by a dot in Fig. |2 These revivals (the first one in particular) can be important 
for experimental applications: in principle they would enable one to obtain a perfect NOT 
transformation, without reaching the adiabatic regime. It is important to notice that this 
result does not depend on the initial state of the system but is a feature of the chosen path 
(see [iJI for details). 

Finally, we emphasize that similar features (and in particular the presence of the revivals 
in the non-adiabatic regime) hold for a large class of gates. For transformations consisting in a 
loop that starts at the pole and is composed of three geodesies given by two arcs of meridians 
and an arc of the equator, enclosing a solid angle to = 7r/2n (n G N*), there is a straightforward 
generalization of Eq. 1)3. 7|1 : 

(2n + l)7r 



r^n) = 2nl7 V16fc%^-1. (3.9) 

This expression is valid provided that the loop is covered at a constant angular speed [see Fig. 
n^b)]: t?arci = ¥'arc2 = '!?arc3 = const. Reversing the orientation of the loop leads to identical 
results. 

IV. MASTER EQUATION FOR A TIME-DEPENDENT HAMILTONIAN 

The interaction between a system and the environment is usually analyzed in terms of a 
master equation. In the standard approach to this problem one assumes that the Hamiltonian 



of the system is time independent (see for instance |2j|). For time dependent Hamiltonians a 



slightly different approach is needed 
time dependent system Liouvillian 



23]. We consider a general Liouville operator with a 



C{t) = £o(t) + XCsB = Csit) 01 + 1^ Cb + XCsb, (4.1) 

where A is the dimensionless coupling constant representing the strength of the noise. The 
evolution of density operator git), describing the system and the environment, is governed by 



the von Neumann-Liouville equation g{t) = C{t) Q{t). We assume that there are no initial 
correlations between system and bath (i.e. the initial state is factorized) and that the bath is 
in equilibrium (e.g. in a thermal state). The main hypothesis in the derivation of a master 
equation is that the typical timescale of the evolution is much slower than the timescales 
characterizing the bath. We shall also assume that the timescale related to the rate of change 
of the system Hamiltonian is the slowest timescale of our problem: this is clearly related to 
the adiabaticity of the evolution. In other words, compared to the bath correlation time, the 
evolution of £5 is always "adiabatic." This is assured by the condition 

TcA « 1, (4.2) 

where Tc is the correlation time of the bath and the energy gap A = min |e„(t) — em(t)| char- 
acterizes the rate of change of £5. Under these conditions one gets 

ait) = [Csit) + X^Tit)]ait), (4.3) 

where a{t) = Tib {Qi^)} is the system density matrix and T{t) a time dependent dissipation 
superoperator. Equation (|4.3|) is the same master equation one would obtain by consider- 
ing Cs{t) "frozen" at time t and evaluating the decay rates and the frequency shifts at the 
instantaneous eigenfrequencies u;{t) = emit) — e„(t) of the system Liouvillian. 

We consider now the physical system described in Sec. lIIII For simplicity let the environment 
affect only the transitions between levels |0) and |e); this is enough for our purposes. The total 
Hamiltonian is Hxit) = H{t) + Hb + XHsb, where H{t) is the system Hamiltonian. The 
bath is bosonic, Hb = "^j^cuko-k o,k with w^ the frequency of the A:-th mode. The interaction 
Hamiltonian is Hsb = J2k1k{\0){e\ + |e)(0|) (a^^ + a^), where 7^ is the coupling constant 
between the system and the k-th mode of the bath. By using Eq. 1)3. 2|) and the form of 
the interaction Hamiltonian, we can obtain time dependent Lindblad operators describing the 
transitions caused by the environment. In the interaction picture generated by the operator R 
defined in 1)2. Ijl . the density operator aji{t) = i?V(t)i? satisfies the following master equation: 

anit) = -i[Hs{0),aRit)] - i[D{t, 0), aR{t)] + X^T{t)aRit), (4.4) 



where the Lamb shifts and the decay rates can be evaluated by standard formulas 
when one introduces the appropriate thermal spectral densities. 
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FIG. 3: Mean fidelity (F) versus cyclic time fir. The dissipation constant A^ increases from top to 
bottom: A^ = (noiseless case), 0.005, 0.01, 0.02, 0.03, 0.04 and 0.05. 

V. FIDELITY AND BEHAVIOR OF THE OPTIMAL WORKING POINT 



Equation (|4.4)) was numerically integrated along the loop in Fig. WO^) when the three arcs 
are covered at a constant angular speed. The values of the Lamb shifts and decay rates were 
assigned, somewhat arbitrarily, for illustrative purposes. They correspond to a bath at very 
high temperature. 

The behavior of the average fidelity ()3.4|) . with a{0,T, •&, if) numerically obtained from 1)4. 4() . 
is shown in Fig. |3| from top to bottom, the dissipation constant increases from A^ = to 0.05. 
In the noiseless case (upmost line) the fidelity tends to 1 when fir — > oo (adiabatic limit). This 
asymptotic value is not reached monotonically: there are some oscillations, with maxima at 
i*" = 1 in the noiseless case. This is the case discussed in Section HTll the NOT transformation 
is perfect, even though one is far from the adiabatic regime, at the time values given by (|3.7j) . 

Clearly, in the presence of noise, the fidelity decreases as the time needed for the transforma- 
tion increases. This can make it difficult to obtain a pure geometrical transformation (because 
of the necessary adiabatic condition). Therefore, it appears convenient to take advantage of 
the presence of the peaks. As a matter of fact, the fidelity decrease due to the noise is very 
small in the non adiabatic regime and one can think of realizing the NOT gate by fine tuning 
the total operation time. As the best performance is obtained for the first peak of the fidelity 
{optimal operation point), we shall focus on the A-dependence of the coordinates of the first 
significant maximum, F* and r*, and their deviation from the noiseless values Fq and rj* (see 




0.0002 0.0004 0.0006 0.0008 0.001 



0.01 0.02 0.03 0.04 0.05 0.06 

A2 



FIG. 4: Maximum value of the mean fidelity F* vs noise (coupling to the bath) A^. (a) Small coupling: 
the fit yields F* = 1 - 6.34A^ (b) Larger coupling: the fit yields F* = 1 - 6MX'^ + 29.93A'^; linear fit 
as in (a). The error bars are always smaller than the size of the points. 



Fig. [2}. It is important to stress that in the non adiabatic regime the gate is no longer pm'ely 
geometrical; in principle it would be possible to extract the geometric contribution, but one 
would not gain any additional information, useful for experimental purposes. 

A critical issue is the total amount of noise. In the simulations in Fig. |31 we considered 
a noise strength A^ ranging from 0.005 to 0.05. However, a realistic physical estimate, using 
thermal spectral densities, would yield a noise level below 0.005 [lj|. In this regime the fidelity 
at the optimal point reaches values greater than 0.9. From this result it is clear that we can 
exploit the optimal times for realizing the NOT transformation with a relatively high fidelity 
even in absence of additional control. 

It is important to understand how the optimal time and the corresponding value of fidelity 
change by increasing the strength of the noise. This should yield information about the ro- 
bustness of holonomic quantum computation against the detrimental effects of noise. Figure 
@fa) shows the behavior of F* for small noise: the points are the result of a numerical analysis 
and the continuous line is the fit 



1 - F2X', 



(5.1) 



yielding F2 = 6.34. The agrement is excellent and enables one to conclude that fidelity de- 
creases linearly with A^ for A^ < 10~^. The behavior of F* for larger values of A^ is displayed 
inFig.Etb). The fit is 



F* = 1- F2A2 + F4A^ 



(5.2) 



with F4 = 29.93. Observe that A^ ~ 6 x 10 ^ is a very large (unphysical) value. The conclusion 
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FIG. 5: Optimal time r* vs noise (coupling to the bath) A^. The values of r* are more affected by 
error than those of F* . (a) SmaU couphng: the fit yields fir* = firj* — 59.40A^, with firj* = 18.25, the 
theoretical value (ESJ. (b) Larger couphng: the fit yields nr* = f^rj* - 59.40A2 + ggO.eSA"* - 7655. OSA^; 
linear fit as in (a). 



that fidelity decreases as A^ (for small A) is to be expected from a perturbation expansion of 
the master equation ()4.4|) . 

Let us now analyze the behavior of the optimal time r* vs A^. Figure [H^a) displays the 
otpimal time vs A^ in the small coupling case. The values of r* are more affected by errors due 
to the approximations introduced by the numerical calculations. The optimal time decreases 
by increasing A^. The fit yields 



rsA^, 



(5.3) 



with T2 = 59.40fi^"'^, the value r* = 18.25fi^"'^ being obtained analytically from Eq. H3.8() . 
Again, the linear fit is in good agreement with the data and enables one to conclude that the 
optimal time decreases linearly with A^ for A^ < 10^^. The behavior of r* for larger values of 
A^ is displayed in Fig. [5I^b). In this case the fit is 



r* = ri* - raA^ + r4A^ - rgA^ 



(5.4) 



with T4 = 990.65fi"^ and tq = 7655. 95fi"^ 

The behavior of the point (F*,t*) is shown in Fig. El both for small (a) and larger (b) 
values of A^. From Eqs. ()5.1|) and (|5.3|) one obtains 



F* = 1- -^n{T^ -T*) = l + 0.11fi(r* - ri ' 
fir2 



(5.5) 



This is the small coupling situation displayed in Fig. El^a). Notice that the mean fidelity 
increases linearly with r* in this regime. In the presence of noise, optimal quantum gates are 



11 



p* 



1 



0.999 
0.998 
0.997 
0.996 
0.995 
0.994 
0.993 




18.19 18.2 18.21 18.22 18.23 18.24 18.25 
flT* 



FIG. 6: Maximum value of the mean fidelity vs the optimal working time t*. (a) Small coupling: the 
fit yields F* ~ 1 + 0.11il(r* — rj*). (b) Larger coupling; the errors are included in the size of the points 
and the linear fit is as in (a). 
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FIG. 7: Mean fidelity (F) versus cyclic time fir: a closer look at the evolution of the optimal point in 
Fig. 121 The dissipation constant A^ increases from top to bottom: A^ = (noiseless case), 0.005, 0.01, 
0.02, 0.03, 0.04 and 0.05. The dotted line is the fit described in the text. 



less precise and slightly faster than the ideal ones. A further analysis of this dependence for 
larger values of the noise [using Eqs. (|5.2() and H5.4() ] yields an involved algebraic expression 
that includes higher order corrections. The behavior is shown in Fig. El^b). The evolution of 
the optimal working point is summarized in Fig. [3 
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FIG. 8: Robustness parameter R vs X^. 
VI. ROBUSTNESS AND CONCLUDING REMARKS 

In order to shed some light on the robustness of the quantum gate in the neighborhood of the 
optimal working point it is useful to compare the optimal fidelity F* with the fidelity obtained 
in the adiabatic limit -Fadiab- Let us observe that in general, the analysis of decoherence in 
geometric computation raises a critical issue in connection with adiabatic evolutions, that 
cannot be too slow, as decoherence would eventually destroy any interference. This inevitably 
introduces an element of arbitrariness in the following definitions. 

We shall evaluate the performance of the optimal (non adiabatic) quantum gate, as com- 
pared to its adiabatic limit, throught a "robustness" parameter: 

R=- ^^. (6.1) 

By glancing at Figs. |21 and 01 it is apparent that the adiabatic limit is practically attained 
already at the third peak, namely for r = r| in Eq. (|3.7j) . We therefore take 

The dependence of R on A^ is displayed in Fig. |H1 Clearly, the (relative) robustness of the 
optimal gate is larger for larger noise levels. We notice the presence of a linear regime for small 
coupling. 

Although we focused our attention on the particular physical system shown in Fig. ^a), 
our conclusions are rather general and are valid for other physically relevant situations. There 
exist some values of the duration of the evolution for which the fidelity is 1, even though one 
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is far from the adiabatic regime. As already emphasized at the end of Sec. II I II these results 
can be extended to more general loops, yielding optimal times like in Eq. (jSIHI)- 

The presence of these optimal peaks is important for experimental applications: if the total 
operation time can be tuned to the first peak, one can realize a transformation that is the most 
robust against noise. Moreover the maximum value of the fidelity exhibits different regimes: for 
small coupling it decreases linearly with A'^ . In general, the gate can be considered robust when 
compared to the standard adiabatic one. This can be of interest for experimental applications, 
if one aims at introducing further control. The case of two-qubit gates is not trivial and is at 
present under investigation. 
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